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1. Introduction 
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Figure 1: Three examples of the modern AW potential in 1 Sq (spin singlet and 5-wave) channel: Bonn[||], 
Reid93[0 and AV18[§. 

In 1935, in order to explain the origin of the nuclear force which binds protons and neutrons 
(nucleons) inside nuclei, Yukawa introduced virtual particles, % mesons, an exchange of which be- 
tween nucleons produces the famous Yukawa potential [[T]]. Since then, both theoretically and exper- 
imentally, enormous efforts have been devoted to understand the nucleon-nucleon (AW) potential, 
recent examples of which are displayed in Fig.|l|. These modern AW potentials are characterized 
as follows[Q, ||]. At long distances (r > 2 fm) there exists weak attraction generated by the one 
pion exchange potential(OPEP), and contributions from the exchange of multi-pions and/or heavy 
mesons such as p makes an overall attraction a little stronger at medium distances ( 1 fm < r < 2 
fm). At short distances (r < 1 fm), on the other hand, attraction turns into repulsion, and it becomes 
stronger as r becomes smaller, forming the strong repulsive core|Q]. Although the repulsive core is 
essential not only for describing the AW scattering data but also for the stability of atomic nuclei, 
its origin remains one of the most fundamental problems in nuclear physics for a long time[||]. 

In a recent paper[]|], using lattice QCD simulations, three of the present authors have calculated 
the AW potential, which possesses the above three features of the modern AW potentials, as shown 
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Figure 2: The central (effective central) AW potential (MeV) as a function of r (fm) for the singlet (triplet) 
at m K ~ 530 MeV in quenched QCD[g]. 
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in Fig.^| . This result has received general recognition[[T(]] . 

The above potentials have been extracted from the Schrodinger equation as 

V E (x)(p E (x) = (e + ^-\ <p E (x) (1.1) 

with the reduced mass m = mjy/2, using the equal-time Bethe-Salpeter wave function (p E (x), de- 
fined by 

(p E (x) = (0\N(x,0)N(0,0)\NN;E), (1.2) 

where \NN;E) is an eigen-state of two nucleons with energy E and N(x,t) is an interpolating 
operator for the nucleon. The potentials in Fig.|| are obtained at E ~ 0. From this definition it is 
clear that the potential V E (x) may depend on the value of energy E and/or the choice of the operator 
N(x,t). In this talk, we focus on the energy dependence of the potential V E (x). In Sect.||, V E (x) 



from an integrable model in 2 dimensions is considered! H ]• The AW potential calculated at E ^ 



in quenched QCD is presented in Sect.[3| Our discussions are given in Sect.||. 

2. Potentials from an integrable model 

In this section we consider the Ising field theory in 2 dimensions, where the one-particle state 
of mass M and the rapidity 6, denoted by \6), has momentum p = M(cosh0,sinh0) with state 
normalization 

(6'\6) = 4n8(6-d'). (2.1) 
The Bethe-Salpeter wave function is defined by 

W(r,0) = /(O|a(^,O)a(O,O)|0,-0) in (2.2) 
where \6, — 6} m is the 2-particle in-state and r = Mx. The spin field a(x) is normalized as 

{0\o(x)\d)=e-*\ (2.3) 



The explicit form of this wave function has been calculated by Fonseca and Zamolodchiko v Q 1 2|] as 

e X(r)/2 



cosh0 

where <I>±, q> and % satisfy 



+ (r, df cosh -0j- 3>-(r, d) 2 cosh + d 



(2.4) 



4>' ± (r,0) = ^sinh(<p(r)±e)O T (r,0), (2.5) 

\Wir)]' = lsmh(2<p(r)), (2.6) 

- r [rX'{r)]' = i[l-cosh(2<p(r))]. (2.7) 
In the limit that r — > 0, the wave function has the expansion 

*F(r,e) ~Cr 3/4 sinh(0) + O(r 7/4 ), (2.8) 
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Figure 3: The Ising potential (multiplied by r 2 ) r 2 V e (r) for 6 = 1 .0 (dotted), 9 = 0.6 (dot-dashed), = 0.3 
(dashed) and 9 = (solid). 



which is expected from the operator product expansion (OPE), 

ff(jc,0)ff(0,0) ~ G(r)l + cr 3/4 £ (0) + • • • . 
where S{x) is the mass operator of dimension 1. 



(2.9) 



We can solve the coupled equations (2.5), (2.6) and (2.7) for <I>±, ip,^ numerically with their 
boundary conditions at r = 0[11|. From the wave function, a rapidity-dependent potential can be 
obtained by 

x ¥"(r,d) + sinh 2 d x ¥(r,d) 



Ve(r) 



(2-10) 



As r 



(2.11) 



0, however, from fl2.8|), the potential Vq becomes rapidity independent: 

where not only the power of r is universally -2 but also the overall coefficient is determined as 3/4 x 
(3/4— 1) from the r 3 / 4 behaviour of the wave function. In Fig.[3], r 2 Vg(r), the potential multiplied 
by r , is plotted as a function of r for several values of 6. We observe that an energy (rapidity )- 
dependence of potentials is small at 6 < 0.6 1 . In particular, potentials are almost identical between 
6=0 and 6 = 0.3. The energy dependence of the Ising potential seems weak at low energy. 
Although the physics in the Ising model is vastly different from QCD, we hope that a similar 
property holds for the AW potential. In the next section we investigate an energy dependence of the 
AW potentials in quenched QCD. 



3. Nucleon-nucleon potentials at non-zero energy in quenched QCD 



We follow the strategy in Ref.[|9|, |13J to define the wave function and to calculate the potential 
through it. Let us explain our set-up of numerical simulations. Gauge configurations are gener- 
ated in quenched QCD on a 32 3 x 48 lattice with the plaquette gauge action at j8 = 5.7, which 

'Note that the singularity of the potential for 6 = 1.0 is caused by the vanishing of the corresponding wave function 
at this point. 
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S n NN Wave Function 
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Figure 4: The AW wave function with APBC (red bars) at t — 11a, together with the fit by the Green's 
function at large distances (green crosses). 
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Figure 5: Left: The central AW potentials for the 1 Sq state with APBC (red bars) and PBC (blue crosses) in 
quenched QCD at t = 9a. Right: Its zoom-in. 



corresponds to a ~ 0.137 fm. We employ the Wilson quark action with anti-periodic boundary 
condition (APBC) in space at the hopping parameter K = 0.1665, corresponding to m n ~ 530 MeV 

71 

and niN — 1330 MeV. The minimum momentum is given by p m j n = (1, 1,1), which leads to 

32a 

k 2 

Pmin — 240 MeV and E = — ~ 50 MeV, where E is the non-relativistic energy in the center of 

m N 

mass system. 2000 configurations are accumulated to obtain our result. 

We first determine a value of k 2 , by fitting the wave function at large distance (13a < |x| < 16a) 
with the Green's function of the Helmholtz equation on an L 3 box, given by 



n x ,n y ,n z GZ ,(3.1) 



1 e i(2n/L)a-x f / 1 1 1 

G(x; * } = i?]* (2x/L)he-* ' r= 1 { n *+r n y + 2> n * + 2 

as plotted in Fig.^. The fit gives k 2 a 2 = 0.030(4) with # 2 /dof ~ 1.8 at t = 11a , which corresponds 
to E ~ 50 MeV in physical unit. 

In Fig.|, the central NN potential for the 1 So state with APBC ( E ~ 50 MeV) is plotted as a 
function of r at t = 9a, together with the one with PBC ( E ~ ). Fluctuations of data with APBC 
at large distances ( r > 1 .5 fm ) are mainly caused by contaminations from excited states, together 
with statistical noises. Data at larger t are needed to reduce such contaminations from excited 
states, though statistical errors also become larger. A non-trivial part of potential at r < 1.5 fm, on 
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the other hand, are less affected by such contaminations. As seen from Fig. the AW potentials 
are almost identical between E ~ and E ~ 50 MeV. 

4. Discussion 

As discussed in the introduction, the potential defined from the Bethe-Salpeter wave function 
depends on the energy: 

V E (x)<p E (x) = (e + ^-\ (p E {x). (4.1) 



In Q13 , 14], it is shown that the energy-dependent potential can be converted to the energy-independent 



but non-local potential as 

J d 3 yU(x,y)(p E (y) = (e + ^) <p E (x) = V E (x) <p E (x) . (4.2) 

We then apply the derivative expansion to this non-local potential [|T3|] as 

U(x,y) = V(x,V)8(x-y) (4.3) 
V(x,V) =Vo{r)+V a {r){a i -a 2 ) + V T {r)Sn + 0{W) (4.4) 

where r = |x|, 0\ t 2 represents the spin of nucleons, and 

3 

5i2 = -^(o l -x)(o 2 -x)-(o l -o 2 ) (4.5) 

is the tensor operator. Our result in the previous section indicates that non-locality is very weak. 

The analysis for the potentials in the Ising field theory in 2 dimensions suggests an interesting 
possibility that the universality of potentials at short distance can be understood from a point of 
view of the operator product expansion (OPE). If this is the case, the origin of the repulsive core 
might be explained by the OPE. We are currently working on this problem. 

Before closing this talk, we consider an alternative possibility to construct the energy-independent 
local potential [jl4"l]. The inverse scattering theory suggests that there exists an unique energy inde- 
pendent potential, which gives the correct phase shift at all energies. Here we propose a new 
method to construct the energy-independent local potential from V E . For simplicity , the 1 dimen- 
sional case is considered. Suppose that <l>(x) = A E (x)(p E (x) satisfies the Schrodinger equation with 
the energy-independent local potential V{x), 

d 2 \ 

— + V{x)\ (A E (x) <Pe {x))=E (A £ (x) <p E (jc) ) , (4.6) 

we obtain the following differential equation, 

V(x)A E (x) = V E {x)A E {x) + A' E (x) + 2A^(x)(log (p E {x))' . (4.7) 

(Here we set 2m = 1 .) If V(x) is given, A E (x) can be easily obtained from this equation. We first 
consider a finite box with size L, which allows only discrete momenta, k n ~ 27in/L, n = 0, 1,2, • • • . 
Once E n = k„ is given, q> E (x) becomes zero at n + 1 points £l n = {xq,x\ , • • • ,x n }. We then have 

= K E (xi)A E (xi) +2A , E (x i )<p' E (x i ) (4.8) 
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for Xi € Q. n , where Ke(x) = (E + d 2 /dx 2 )q>E(x). Since Q. n becomes dense in [0,L] in the n — > oo 
limit, V(jc) can be constructed as 

V (x) = lim { V E (jc) - 2X £ (jc) (log (p E (x) )' - X' E (x) + X E (x) 2 } (4.9) 

K E (xi) 



where X E (x) is an interpolations of Xe(xi) = - — — — r- with x, G O n . If the limit ( p^9| ) exists, the 

energy-independent local potential V (jc) can be obtained. In the 3 dimensional case, we first intro- 
duce the polar coordinate, and then apply the above procedure in 1 dimension to the radial variable 
r with the fixed angular momentum /. 
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